Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,449,937 G→T 6.4% R99R (CGC→CGA rplV ← 50S ribosomal subunit protein L22

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,449,9370GT 6.4% 342.5 / 6.9 93R99R (CGC→CGArplV50S ribosomal subunit protein L22
Reads supporting (aligned to +/- strand):  ref base G (41/46);  new base T (3/3);  total (44/49)
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.30e-01

AATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCA  >  NC_000913/3449798‑3450075
                                                                                                                                           |                                                                                                                                          
aaTTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATAGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGAGATGTGGCTGGTg                                                                                                                                            >  2:785131/1‑140 (MQ=255)
  ttCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACCGTGATGTGGCTGGTgcg                                                                                                                                          <  2:593737/140‑1 (MQ=255)
    cTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCt                                                                                                                                        <  2:736632/140‑1 (MQ=255)
      ttGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTc                                                                                                                                      <  2:803554/140‑1 (MQ=255)
      ttGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTc                                                                                                                                      <  2:746158/140‑1 (MQ=255)
       tGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCa                                                                                                                                     >  1:358725/1‑140 (MQ=255)
        ggTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAg                                                                                                                                    >  1:1732872/1‑140 (MQ=255)
           gtTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGAt                                                                                                                                 <  1:107221/140‑1 (MQ=255)
           gtTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGAt                                                                                                                                 <  1:107228/140‑1 (MQ=255)
            ttCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCCGAGTCTCAGCGAGCGGACACAACCACAGTGATGTGGCGGGTGCGCTTCAGGATg                                                                                                                                >  1:1503527/1‑140 (MQ=255)
            ttCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATg                                                                                                                                >  1:1503537/1‑140 (MQ=255)
             tCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGGTGGTGCGCTTCAGGATGc                                                                                                                               >  2:1389442/1‑140 (MQ=255)
                   aCCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTg                                                                                                                         >  1:555506/1‑140 (MQ=255)
                     cAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCa                                                                                                                       <  1:364350/140‑1 (MQ=255)
                        gTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGa                                                                                                                    >  1:784414/1‑140 (MQ=255)
                        gTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGa                                                                                                                    >  1:785420/1‑140 (MQ=255)
                             gTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCttt                                                                                                               >  1:1706696/1‑140 (MQ=255)
                             gTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCttt                                                                                                               >  1:1705443/1‑140 (MQ=255)
                                 cATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCa                                                                                                           <  2:1379636/140‑1 (MQ=255)
                                  aTGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAc                                                                                                          >  1:503921/1‑140 (MQ=255)
                                  aTGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAc                                                                                                          >  1:220528/1‑140 (MQ=255)
                                  aTGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAc                                                                                                          <  2:718259/140‑1 (MQ=255)
                                   tGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAcg                                                                                                         <  1:1245782/140‑1 (MQ=255)
                                   tGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAcg                                                                                                         <  1:120993/140‑1 (MQ=255)
                                   tGGTTTTACAATACCCAGGCGAATACCATTAGGAGGTACTTTCTGACCCATTGCTAGTCTCCAGTGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAcg                                                                                                         <  1:686592/140‑1 (MQ=255)
                                     gTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAcgcg                                                                                                       <  1:839224/140‑1 (MQ=255)
                                           ccaTACCCCGGCGAATACCACTAGGCTGAACTTTCTGCCCCATTGCCAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGTTGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATaa                                                                                                 <  2:21188/138‑1 (MQ=255)
                                                cccAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATgcgc                                                                                            >  2:558023/1‑140 (MQ=255)
                                                cccAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATgcgc                                                                                            >  1:13502/1‑140 (MQ=255)
                                                cccAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATgcgc                                                                                            >  1:1094629/1‑140 (MQ=255)
                                                cccAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATgcgc                                                                                            >  2:560175/1‑140 (MQ=255)
                                                 ccAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCt                                                                                           <  2:1642409/140‑1 (MQ=255)
                                                    ggCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCa                                                                                        <  2:1566018/140‑1 (MQ=255)
                                                     gCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCAt                                                                                       >  2:1817740/1‑140 (MQ=255)
                                                        aaTACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCt                                                                                    >  1:895028/1‑140 (MQ=255)
                                                        aaTACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCt                                                                                    >  1:114615/1‑140 (MQ=255)
                                                        aaTACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCt                                                                                    >  1:895220/1‑140 (MQ=255)
                                                         aTACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTc                                                                                   <  1:10611/140‑1 (MQ=255)
                                                         aTACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTc                                                                                   >  1:1539945/1‑140 (MQ=255)
                                                          tACCATTAGGATGTCCTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGTTCGGACTCAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCCCGTCCTTTTGCTCGCGGCAGAATGCGCTTCATGcgcg                                                                                  >  2:1683233/1‑140 (MQ=255)
                                                          tACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCg                                                                                  >  2:1405522/1‑140 (MQ=255)
                                                                       tACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTAc                                                                     >  1:1472743/1‑140 (MQ=255)
                                                                       tACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCAGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCATCGTCTAc                                                                     >  2:589397/1‑140 (MQ=255)
                                                                       tACCTTCTTACCCATTGCTAGTCACCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTTGCTGGTGCGCTTAAGGATGCGATCTGCACGACCGTGTGCACGCCGCATAATGCGCTTCATGCTCGGGCCGACGTCTAc                                                                     >  2:589396/1‑140 (MQ=255)
                                                                            tCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGaaaa                                                                >  1:1708188/1‑140 (MQ=255)
                                                                            tCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGaaaa                                                                >  2:503121/1‑140 (MQ=255)
                                                                            tCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGaaaa                                                                >  1:871486/1‑140 (MQ=255)
                                                                            tCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGaaaa                                                                >  1:1708285/1‑140 (MQ=255)
                                                                            tCTGACCCATTGCTACACTCCAGAGTCTCAGCGATCGGACACAACCACAGAGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTATGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGaaaa                                                                >  1:770149/1‑140 (MQ=255)
                                                                               gACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTTCGCTTCAGGATGCGATCTGCACTACCTTTTGCATGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAAttt                                                             >  2:1263678/1‑140 (MQ=255)
                                                                               gACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTTCGCTTCAGGATGCGATCTGCACTACCTTTTGCATGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAAttt                                                             >  2:1263662/1‑140 (MQ=255)
                                                                               gACCCATTGCTAGTCTCCAGAGTCCCAGCGATCGGACACAACCACAGTGATGTGGCTGGTTCGCTTCAGGATGCGATCTGCACTACCTTTTGCATGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAAttt                                                             >  2:1264501/1‑140 (MQ=255)
                                                                                   cATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGt                                                         <  1:1817740/140‑1 (MQ=255)
                                                                                   cATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGt                                                         <  1:272669/140‑1 (MQ=255)
                                                                                      tGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAAc                                                      >  1:1411564/1‑140 (MQ=255)
                                                                                           gTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCa                                                 >  2:962488/1‑140 (MQ=255)
                                                                                             ctcCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGa                                               <  1:1395457/140‑1 (MQ=255)
                                                                                             ctcCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGa                                               <  1:1010476/140‑1 (MQ=255)
                                                                                             ctcCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGa                                               <  1:661888/140‑1 (MQ=255)
                                                                                                  gagTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTTCGCTTCAGGATGCGATCTGCACTACCTTTTGCATGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTc                                          <  1:1263678/140‑1 (MQ=255)
                                                                                                  gagTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTTCGCTTCAGGATGCGATCTGCACTACCTTTTGCATGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTc                                          <  1:1263662/140‑1 (MQ=255)
                                                                                                  gagTCCCAGCGATCGGACACAACCACAGTGATGTGGCTGGTTCGCTTCAGGATGCGATCTGCACTACCTTTTGCATGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTc                                          <  1:1264501/140‑1 (MQ=255)
                                                                                                   agTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCa                                         <  2:503865/140‑1 (MQ=255)
                                                                                                   agTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCa                                         <  2:936611/140‑1 (MQ=255)
                                                                                                     tctcAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAAt                                       >  1:60921/1‑140 (MQ=255)
                                                                                                       tcAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTACATGCTCGGGCCTTCGTCTACGAAAATTTTCGAAACTGTCAGATCGTCAATGt                                     >  2:649278/1‑140 (MQ=255)
                                                                                                        ccgCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGTCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTc                                    <  2:1381147/138‑1 (MQ=255)
                                                                                                        cAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTc                                    >  2:674774/1‑140 (MQ=255)
                                                                                                          gCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAg                                  <  1:1546938/140‑1 (MQ=255)
                                                                                                          gCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAg                                  <  1:1440941/140‑1 (MQ=255)
                                                                                                          gCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAg                                  <  1:1440969/140‑1 (MQ=255)
                                                                                                             aTCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAgcgc                               <  1:482226/140‑1 (MQ=255)
                                                                                                               cGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCa                             <  1:346768/140‑1 (MQ=255)
                                                                                                               cGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCa                             <  1:756252/140‑1 (MQ=255)
                                                                                                               cGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCa                             <  1:346909/140‑1 (MQ=255)
                                                                                                               cGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGAGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTAAATGTCAGCGCCa                             <  1:1080046/140‑1 (MQ=255)
                                                                                                                   cacaACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGt                         <  1:245415/140‑1 (MQ=255)
                                                                                                                    acaACCACAGTGATGTGGCTGGTGCGCTTCAGGAAGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGtt                        <  2:1660759/140‑1 (MQ=255)
                                                                                                                     caACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTtg                       <  2:1077498/140‑1 (MQ=255)
                                                                                                                       aCCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTtgtg                     >  1:762729/1‑140 (MQ=255)
                                                                                                                         cacaGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGtt                   <  2:762729/140‑1 (MQ=255)
                                                                                                                         cacaGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGtt                   <  1:63849/140‑1 (MQ=255)
                                                                                                                         cacaGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGtt                   <  1:92280/140‑1 (MQ=255)
                                                                                                                         cacaGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGtt                   <  1:1465104/140‑1 (MQ=255)
                                                                                                                            aGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAg                >  1:664098/1‑140 (MQ=255)
                                                                                                                              tGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCg              <  2:789702/140‑1 (MQ=255)
                                                                                                                                aTGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCCATGTCAGCGCCATCGTTGTGTTCAGCGtt            >  2:852735/1‑140 (MQ=255)
                                                                                                                                aTGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGtt            >  2:1050750/1‑140 (MQ=255)
                                                                                                                                  gtgGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAg          <  2:4083/140‑1 (MQ=255)
                                                                                                                                    ggCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCa        <  2:531877/140‑1 (MQ=255)
                                                                                                                                    ggCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCa        <  2:1164686/140‑1 (MQ=255)
                                                                                                                                    ggCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCa        <  2:1378233/140‑1 (MQ=255)
                                                                                                                                    ggCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCa        <  2:1378284/140‑1 (MQ=255)
                                                                                                                                    ggCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCa        <  2:1418795/140‑1 (MQ=255)
                                                                                                                                          tGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCa  >  2:45013/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                          
AATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCA  >  NC_000913/3449798‑3450075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: