Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,530 2 bp→TT coding (114‑115/1260 nt) rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,5300CT97.7% 130.2 / ‑6.3 43I38I (ATC→ATT) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/19);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
*NC_0009133,966,5310CT97.6% 123.8 / ‑4.4 43L39L (CTG→TTG) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/18);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATCCTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGT  >  NC_000913/3966298‑3966772
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                     
cTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑Caa                                                                                                                                                                                                                                      >  1:194036/1‑249 (MQ=255)
    gACCAGGCAGCGAAAAGACATGAGTCGAGGACCGTAAACAGGCAGGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTg                                                                                                                                                                                                                                <  2:505222/251‑1 (MQ=255)
          gCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGGGGCGa                                                                                                                                                                                                                            >  1:387138/1‑249 (MQ=255)
              ggAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGata‑t                                                                                                                                                                                                                     <  1:369473/250‑1 (MQ=255)
                 aaaGACATGAGTCGAGGAACGAAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCtt                                                                                                                                                                                                                  <  2:330889/251‑1 (MQ=255)
                          aGTCGATGACCGTAAACAGGCCTGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGACGAAAATATGGGGCTGGACAACCTGGCTCTTTTGCGTACGCCGGACATTATTTTTTCC‑CTTTTGAA‑GCAGCACG‑CAAAACGTGGCGAAGATA‑GCTTTGGttattg                                                                                                                                                                                                         >  2:76196/1‑246 (MQ=255)
                           gTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGc                                                                                                                                                                                                        >  1:484874/1‑251 (MQ=255)
                                        aaCAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGatat                                                                                                                                                                                           <  1:474545/251‑1 (MQ=255)
                                          cAGGCATGGATGAGCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCGGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACGGGAGATATTg                                                                                                                                                                                         <  1:76196/251‑1 (MQ=255)
                                           aGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGc                                                                                                                                                                                        >  1:257120/1‑251 (MQ=255)
                                                 ggATGATCCTGCCATACCAGTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCTCTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAgga‑tg                                                                                                                                                                                 <  2:460515/251‑1 (MQ=255)
                                                  gcggatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAACGAATACGCCGGTTTCTGAGCTGATCACTCTCTGCGTAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑TTTTTGAC‑GCAGCACT‑CAAAGATTGGCGAAGATA‑TCTTTTGTGTTGGCGTACTGGAGATATTGCAgga‑tgg                                                                                                                                                                                >  2:145861/4‑251 (MQ=255)
                                                     gacccTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGAtt                                                                                                                                                                             <  2:387138/248‑1 (MQ=255)
                                                      atCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACTACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGAttt                                                                                                                                                                            >  2:371868/1‑251 (MQ=255)
                                                           gCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACGCACCCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACGCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGGATGCGTAAGCAGGACATTATTGTTGCC‑AGTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATGTGGtt                                                                                                                                                                        <  2:95724/250‑1 (MQ=255)
                                                                aCCATTCACAGCATTAGGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATAACTCTCGGCGAAAATATGAGGCTGGAAAACCTGGCTCGTATGCGTACGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAAATA‑TCTTTGGTGATGTGGTACTGGAAATATTGACGGA‑TGGATTTGGTTtcctcc                                                                                                                                                                  >  2:95457/1‑251 (MQ=255)
                                                                   aTTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGtt                                                                                                                                                               <  1:10160/251‑1 (MQ=255)
                                                                    ttCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGC‑TGGATTTGGTTTCCTCCGTTc                                                                                                                                                              >  1:492121/1‑251 (MQ=255)
                                                                    ttCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGGGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGttt                                                                                                                                                              >  2:309150/1‑250 (MQ=255)
                                                                          acaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGGGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGc                                                                                                                                                           >  1:430631/1‑248 (MQ=255)
                                                                            aTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAACGAGTGGCGAAAGTA‑TCTTTGGTGATGGCGGACTGGAGGTATTGCAGGG‑TGGGTTTGGTTTCCTCCCTTCCgaggacc                                                                                                                                                      >  2:133033/1‑245 (MQ=255)
                                                                            aTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGc                                                                                                                                                                                                        <  1:80736/202‑1 (MQ=255)
                                                                                     gagaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGGATGCGTAAGCAGGACATTCTTTTTGCC‑ATCCTGAA‑TCAGCACG‑CAAAAAGTGGCGACG‑TACTCTTTGGTGATGTCGTACCGGAGATATTGCGGGG‑TGGGTTTGGTTTCCTCCTTTGCGCCGACaacgcctagc                                                                                                                                             >  2:475655/1‑242 (MQ=255)
                                                                                                       aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGAc                                                                                                                            >  2:364648/1‑250 (MQ=255)
                                                                                                          aaCTCACACCACTATGAATCTTCCCGCATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACCTTATTTTTGCC‑ATTTTGAC‑GCAGCACG‑AAAAGAGTGGCGCAGATA‑TCTTTGGTGATGGCGTACTGGAGAATTTGCAGGA‑GGGATTTGGTTTCCTCCGTTCCGCCGAACGCTCCTCCCTCGCCGGTCCTGTTGACATcc                                                                                                                        >  2:470089/1‑250 (MQ=255)
                                                                                                              cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGTTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGCGTGGCGAAGATA‑TCTTTGGTGATGGCGTTCTGGAGATATTGCAGGA‑TGGATTTGGTTGCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGt                                                                                                                    <  2:198639/251‑1 (MQ=255)
                                                                                                                    acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcccc                                                                                                              <  1:230069/251‑1 (MQ=255)
                                                                                                                       aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCGGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTGGGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCt                                                                                                             <  2:40952/249‑1 (MQ=255)
                                                                                                                             cTTACCGAATTAAAGAATACGCCGGGTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGACGATA‑TCTTTGGTGATGGCGTACTTGAGATATTGCAGGG‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCCCCTACCTCGCCGGTCCTGAGTACATCTCCGTTTCCCCAAGCCaaa                                                                                                      >  2:315133/1‑250 (MQ=255)
                                                                                                                                                  ccGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTACGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGCTA‑TCTTTTGTGATGGCGTACTGGGGATATTGCAGGA‑TGGGTTTGGGTTCCTCCGTTCCGCAGCCAGCTCCTCCCTCGACGGCCCTGATGACATCTACGTTTTCCCTAGCCAAATCCGCCCTTTCacccctccgac                                                                                >  2:458689/1‑241 (MQ=255)
                                                                                                                                                          tGAGCTGATCACTATCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGAC‑AGTTTGAA‑GCAGCACG‑CAACGAGTGGCGAAGATA‑TATTTGGTGATGGCGTACTGGGGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCTCAGAAAGCTCCTACCTCGCCGGTCCTGATGGCCTCCACGTTTACCCTAGCCCAATCCTCCGGGTCAACCGCCGCACTGGGGAt                                                                         >  2:257932/1‑250 (MQ=255)
                                                                                                                                                              gtGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATGAGTTTTGCC‑ATTTTGAA‑TCAGCACT‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGAGGGCGTATTGGATATATTGCAGGA‑TGTATTTGGTTTTCTCTGTTCCGCAGACGTCTCCTACCTCGCCGGTCCTTATGACAGCTACTTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAt                                                                    <  1:487910/250‑1 (MQ=255)
                                                                                                                                                                gATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtct                                                                  >  1:169481/1‑251 (MQ=255)
                                                                                                                                                                    aCTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                              >  2:423095/1‑251 (MQ=255)
                                                                                                                                                                         cGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TTTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTaa                                                           <  1:423628/249‑1 (MQ=255)
                                                                                                                                                                            ggAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAACCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                      <  1:443876/250‑1 (MQ=255)
                                                                                                                                                                               aaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTACAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                                                   >  1:492762/1‑251 (MQ=255)
                                                                                                                                                                                            gAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCCAATCCGCCGTTTCAACCTCCGCACTGGTGATAACATCTCTGGTAAGTTTCGCCCGCCGCAAGAAgg                                      >  2:454574/1‑251 (MQ=255)
                                                                                                                                                                                                cccgggcacGTATGGGTAAGCAGGAAATTATTTTTGCC‑ATTTGGAA‑GAAGCACG‑AAAAGGGTGGCGAAGATA‑TCTTTGGTGATGGCGTCCGGGAGATATTGCAGGA‑TGTATTTGGTTCCCTCCGTTCCCCAGACAGCTCCTCCCTCGCCGGCCCTGATGACATCTCCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGaa                                  <  1:326056/243‑1 (MQ=255)
                                                                                                                                                                                                 acTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATTTTTAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATAGCGGACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAACTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTTGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTggac                                 >  2:332008/2‑248 (MQ=255)
                                                                                                                                                                                                            tGCGTAAGCAGTACATTATTTT‑GCCTATTTTAAA‑GCAGC‑CGTCAAAGTGGGGCGAAGATT‑TCTTTGTTGATGGTTTACGGGAGATATTGCAGTA‑GGGATTTTGTTCCCCCCGCCCCGCAAACCCCCCCTCCCTCGCCGGCCCTGTTCCCTCCCCCGTTCCCCCTAGCCAATCCCGCCGTTTCAACCTCCGCACTGGTCATCCCATCTCTGGTAAGATACGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                       <  1:240926/250‑1 (MQ=255)
                                                                                                                                                                                                                    aagtaCATTATTTT‑GCC‑ATTTTGAAAGAAGCACG‑AAAAGAGTGGCGAAGATT‑TTTTTGGGGTGGGCTTTCGGGAAATATTG‑AGGAATGGTTTTGGTTCCCTCCGTCCCGCAAACACCTCTTCCCCCGCCGTTCCTTATGACATCTACTTTCCCCCAAGCCAAATCCGCCGTTTCACCCCCCGCACTGGTGATACCATCTCTGTTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa              <  1:195962/247‑1 (MQ=255)
                                                                                                                                                                                                                                tttGCC‑ATTTTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGCAAGTTAACGAAGt  >  1:129151/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                     
CTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC‑ATCCTGAA‑GCAGCACG‑CAAAGAGTGGCGAAGATA‑TCTTTGGTGATGGCGTACTGGAGATATTGCAGGA‑TGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGT  >  NC_000913/3966298‑3966772

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑