Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 21,550 A→G L48L (TTA→TTG ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091321,5500AG96.9% 227.1 / ‑3.2 64L48L (TTA→TTGribFbifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (27/35);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCG  >  NC_000913/21415‑21683
                                                                                                                                       |                                                                                                                                     
gATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCgg                                                                                                                                   <  1:849459/140‑1 (MQ=255)
 aTACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCggt                                                                                                                                  >  2:700118/1‑140 (MQ=255)
 aTACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCggt                                                                                                                                  >  2:60106/1‑140 (MQ=255)
   aCGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCggtga                                                                                                                                <  1:700118/140‑1 (MQ=255)
   aCGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCggtga                                                                                                                                <  1:60106/140‑1 (MQ=255)
        gCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCggtgatggtg                                                                                                                           <  1:402339/140‑1 (MQ=255)
           tacataATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGgtgatggtgatg                                                                                                                        <  1:42723/140‑1 (MQ=255)
           tacataATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGgtgatggtgatg                                                                                                                        <  1:680741/140‑1 (MQ=255)
           tacataATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGgtgatggtgatg                                                                                                                        >  2:790353/1‑140 (MQ=255)
                aaTCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCtttt                                                                                                                   <  1:790353/140‑1 (MQ=255)
                     cAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAAc                                                                                                              >  1:127166/1‑140 (MQ=255)
                      aGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCggtgatggtgat                                                                                                                         >  1:891408/1‑128 (MQ=255)
                      aGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCggtgatggtgat                                                                                                                         <  2:891408/128‑1 (MQ=255)
                           ggCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAAc                                                                                                        <  2:538614/140‑1 (MQ=255)
                           ggCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAAc                                                                                                        <  2:741926/140‑1 (MQ=255)
                           ggCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTAAACCTCAAc                                                                                                        <  1:954155/140‑1 (MQ=255)
                             ccccGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCa                                                                                                      >  2:550375/1‑140 (MQ=255)
                              cccGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCAc                                                                                                     <  1:429297/140‑1 (MQ=255)
                                 gCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTgg                                                                                                  >  1:354183/1‑140 (MQ=255)
                                        gggTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGtt                                                                                           <  2:812678/140‑1 (MQ=255)
                                          gtgtgtGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTg                                                                                         <  1:851832/140‑1 (MQ=255)
                                          gtgtgtGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTg                                                                                         <  2:127166/140‑1 (MQ=255)
                                          gggtgtgCGGACTATGGGTCAGTGCGACGGCGGGCAGCGCGGTCAGCGCGCGCTGTGACAGGGCGTGCGGGAAGAAGGGCGCAAGCGCAACGTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTg                                                                                         <  2:316001/138‑1 (MQ=255)
                                             tgtgCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTa                                                                                      <  1:556635/140‑1 (MQ=255)
                                                  tGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGAt                                                                                 >  2:617827/1‑140 (MQ=255)
                                                   gACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGaa                                                                                                <  2:208854/124‑1 (MQ=255)
                                                   gACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGaa                                                                                                >  1:208854/1‑124 (MQ=255)
                                                   gACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATa                                                                                >  1:146384/1‑140 (MQ=255)
                                                                ttCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAg                                                                   >  1:481003/1‑140 (MQ=255)
                                                                 tCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGa                                                                  <  1:754695/140‑1 (MQ=255)
                                                                  cGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGAc                                                                 >  2:253624/1‑140 (MQ=255)
                                                                     cGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGa                                                              >  1:107306/1‑140 (MQ=255)
                                                                            cATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCt                                                       <  1:617827/140‑1 (MQ=255)
                                                                                    tCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGaaa                                               >  1:586872/1‑140 (MQ=255)
                                                                                    tCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGaaa                                               <  2:481003/140‑1 (MQ=255)
                                                                                           gcgcTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCg                                        <  1:253624/140‑1 (MQ=255)
                                                                                           gcgcTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCg                                        <  2:586872/140‑1 (MQ=255)
                                                                                           gcgcTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCg                                        <  1:21624/140‑1 (MQ=255)
                                                                                           gcgcTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCg                                        <  1:82672/140‑1 (MQ=255)
                                                                                             gcTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGtt                                      <  2:107306/140‑1 (MQ=255)
                                                                                              cTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTa                                     <  1:550375/140‑1 (MQ=255)
                                                                                                gTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTAcc                                   <  2:945195/140‑1 (MQ=255)
                                                                                                     aGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCa                              <  2:949114/140‑1 (MQ=255)
                                                                                                     aGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCa                              <  2:953489/140‑1 (MQ=255)
                                                                                                       ggCTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCaga                            <  2:146384/140‑1 (MQ=255)
                                                                                                         cTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCa                              >  2:260029/1‑136 (MQ=255)
                                                                                                         cTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCa                              <  1:260029/136‑1 (MQ=255)
                                                                                                         cTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAgt                          >  1:640550/1‑140 (MQ=255)
                                                                                                         cTTGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAgt                          >  2:834097/1‑140 (MQ=255)
                                                                                                           tGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGtt                                      <  2:819075/126‑1 (MQ=255)
                                                                                                           tGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGtt                                      >  1:819075/1‑126 (MQ=255)
                                                                                                           tGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTa                                     >  2:426786/1‑127 (MQ=255)
                                                                                                           tGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTa                                     <  1:426786/127‑1 (MQ=255)
                                                                                                           tGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAgtgt                        >  2:998238/1‑140 (MQ=255)
                                                                                                           tGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAgtgt                        >  1:535632/1‑140 (MQ=255)
                                                                                                           tGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAgt                          >  2:790572/1‑138 (MQ=255)
                                                                                                           tGCAGGAAGAAGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAgt                          <  1:790572/138‑1 (MQ=255)
                                                                                                                gaagaaGGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGt                   <  1:998238/140‑1 (MQ=255)
                                                                                                                  agaagGGCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTg                 >  2:515407/1‑140 (MQ=255)
                                                                                                                      gggCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTa             >  1:810026/1‑140 (MQ=255)
                                                                                                                      gggCGCAAGCGCAACTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTa             >  1:1087768/1‑140 (MQ=255)
                                                                                                                        gcgcaagcgcaaCTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACg           >  2:521774/1‑140 (MQ=255)
                                                                                                                                 caaCTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTgtgcg  >  2:115730/1‑140 (MQ=255)
                                                                                                                                 caaCTTGCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTgtgcg  <  1:834097/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCG  >  NC_000913/21415‑21683

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