Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,029,739 (A)8→7 coding (231/1191 nt) rlmI ← 23S rRNA m(5)C1962 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,029,7320A.97.1% 299.3 / 0.4 68coding (238/1191 nt)rlmI23S rRNA m(5)C1962 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (33/33);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACC  >  NC_000913/1029583‑1029876
                                                                                                                                                     |                                                                                                                                                
cagcagTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGaa                                                                                                                                                 <  1:396048/151‑1 (MQ=255)
 agcagTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGaaa                                                                                                                                                <  2:18023/151‑1 (MQ=255)
 agcagTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGaa                                                                                                                                                 >  2:810439/1‑150 (MQ=255)
     gTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGaaaaaa                                                                                                                                             >  1:683861/1‑150 (MQ=255)
     gTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGaaaaaa                                                                                                                                             >  2:843618/1‑150 (MQ=255)
        gcagcaCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGaaaaaaagca                                                                                                                                         <  2:64571/151‑4 (MQ=255)
          agcaCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAAt                                                                                                                                      <  2:68839/151‑1 (MQ=255)
             aCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTc                                                                                                                                   >  1:913741/1‑151 (MQ=255)
              ccAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCg                                                                                                                                  >  1:610382/1‑151 (MQ=255)
                 gAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGac                                                                                                                                >  2:290656/1‑149 (MQ=255)
                  aaaaTTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAg                                                                                                                              <  2:360793/151‑1 (MQ=255)
                     aTTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGGCGATAGa                                                                                                                             >  2:347012/1‑149 (MQ=255)
                      ttACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTc                                                                                                                          >  1:70401/1‑151 (MQ=255)
                        aCCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCa                                                                                                                         <  2:538088/150‑1 (MQ=255)
                          cGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGAATCAg                                                                                                                        >  1:989760/1‑149 (MQ=255)
                           gAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACg                                                                                                                     >  1:367581/1‑151 (MQ=255)
                            aaaCGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACgg                                                                                                                    <  2:82407/151‑1 (MQ=255)
                            aaaCGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACgg                                                                                                                    >  1:674349/1‑151 (MQ=255)
                             aaCGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACggg                                                                                                                   <  2:588309/151‑1 (MQ=255)
                              aCGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGt                                                                                                                  <  2:969628/151‑1 (MQ=255)
                                   cgatAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCaaa                                                                                                              >  2:207192/1‑150 (MQ=255)
                                     atAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGt                                                                                                           >  1:675472/1‑151 (MQ=255)
                                      tAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGGCAAACGTc                                                                                                          >  1:1049140/1‑151 (MQ=255)
                                       aGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTcc                                                                                                         <  2:910093/151‑1 (MQ=255)
                                           aTACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGAc                                                                                                     <  1:1018330/151‑1 (MQ=255)
                                             aCCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGAcgc                                                                                                   >  2:509337/1‑151 (MQ=255)
                                               ccGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGAcgcgc                                                                                                 <  1:329039/151‑1 (MQ=255)
                                                  gCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGcc                                                                                              >  2:31843/1‑151 (MQ=255)
                                                        cATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGAAATGTCGATAGACTCCGACGGGGCAAACGGCCAGACGCGCGCCCGGGtt                                                                                        >  1:67051/1‑151 (MQ=255)
                                                         aTCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGAttt                                                                                       >  2:500081/1‑151 (MQ=255)
                                                           cAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTg                                                                                      <  2:737715/150‑1 (MQ=255)
                                                                   ccGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCt                                                                               <  1:289460/149‑1 (MQ=255)
                                                                    cGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTgg                                                                             >  2:58633/1‑150 (MQ=255)
                                                                           tAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATaa                                                                      <  2:631973/150‑1 (MQ=255)
                                                                           tAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATaa                                                                      <  2:1049887/150‑1 (MQ=255)
                                                                            aaaCGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAgc                                                                    <  1:612455/151‑1 (MQ=255)
                                                                            aaaCGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAgc                                                                    <  2:923096/151‑1 (MQ=255)
                                                                                   aGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCcgcg                                                              <  2:398144/150‑1 (MQ=255)
                                                                                   aGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGt                                                             <  2:465883/151‑1 (MQ=255)
                                                                                      tGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCt                                                          >  1:132577/1‑151 (MQ=255)
                                                                                       gTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCt                                                          >  2:409941/1‑150 (MQ=255)
                                                                                       gTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCt                                                          >  2:188937/1‑150 (MQ=255)
                                                                                          gAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGTGAATAAGCGCCGCGTGCTAAcc                                                      >  2:990365/1‑151 (MQ=255)
                                                                                           aGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAAcc                                                      >  1:697737/1‑150 (MQ=255)
                                                                                             gCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCAtt                                                   <  1:1057971/151‑1 (MQ=255)
                                                                                              ccATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCAttt                                                  <  1:1048963/151‑1 (MQ=255)
                                                                                              ccATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCAttt                                                  <  2:1022702/151‑1 (MQ=255)
                                                                                               cATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCAtttt                                                 <  1:922854/151‑1 (MQ=255)
                                                                                                aTCTTTTTGAGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGGAAAAAAGCAATGTCGATAGACTCAGACGGGGCAAACGTCCAGAAGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCAtttt                                                 >  1:648938/1‑150 (MQ=255)
                                                                                                 tCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTcc                                               >  1:303111/1‑151 (MQ=255)
                                                                                                    ttttGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTccc                                              <  1:842746/149‑1 (MQ=255)
                                                                                                    ttttGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTcc                                               <  2:774150/148‑1 (MQ=255)
                                                                                                    ttttGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtg                                            >  1:200695/1‑151 (MQ=255)
                                                                                                         cgcCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTgatgat                                       <  1:33581/151‑1 (MQ=255)
                                                                                                            cagccagTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATcaa                                    <  1:590107/151‑1 (MQ=255)
                                                                                                              gccagTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCTGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATcaaca                                  <  1:213799/151‑1 (MQ=255)
                                                                                                                cagTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATcaacaa                                 >  2:57974/1‑150 (MQ=255)
                                                                                                                    cACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtgat                                          <  1:614116/137‑1 (MQ=255)
                                                                                                                    cACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGAATCAGAAGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtgat                                          >  2:614219/1‑137 (MQ=255)
                                                                                                                      cGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGa                           >  2:120442/1‑150 (MQ=255)
                                                                                                                         cATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGTTCAACCGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATATTCCCTGATGATCAACAATATCGATgg                        <  1:766365/150‑1 (MQ=255)
                                                                                                                             tttGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTc                    >  2:1032950/1‑150 (MQ=255)
                                                                                                                              ttGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCa                   <  2:26821/150‑1 (MQ=255)
                                                                                                                                gtgCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCAcc                 <  2:1030270/150‑1 (MQ=255)
                                                                                                                                 tgCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCg                >  2:1062699/1‑150 (MQ=255)
                                                                                                                                  gCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCg                >  2:496574/1‑149 (MQ=255)
                                                                                                                                     tgttgCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAggct           >  2:646106/1‑151 (MQ=255)
                                                                                                                                     tgttgCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAggc            <  1:843473/150‑1 (MQ=255)
                                                                                                                                         gCAAATGGCGGGAAAAAAAACCATGTCGATAGCCTCCGACCGGGCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGGTCAAAAAAATCGATGGTTTCACCGAggctggct       >  1:361349/1‑151 (MQ=255)
                                                                                                                                          cAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCtt      >  2:303684/1‑151 (MQ=255)
                                                                                                                                           aaaCGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCttt     <  1:600106/151‑1 (MQ=255)
                                                                                                                                             aCGGCGGGAAAAAAAACAATGTCGATAGACTCCGACGGGGCAAACGTCCCGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGGGCTAACCATTTTTCCTGATGATCCACAATATCGATGGGTTCACCGAGGCTGGCttt     >  1:281169/1‑149 (MQ=255)
                                                                                                                                             aCGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTAc   <  2:708460/151‑1 (MQ=255)
                                                                                                                                              cggcggGAAAAAAAGCAATGTCGATAGAATCAGACGGGTCAAAAGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTAcc  >  1:954113/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACC  >  NC_000913/1029583‑1029876

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