Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,638,124 T→C V161V (GTT→GTC pitA → metal phosphate:H(+) symporter PitA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,638,1240TC98.6% 240.5 / ‑4.7 70V161V (GTT→GTCpitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (38/31);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGC  >  NC_000913/3637988‑3638260
                                                                                                                                        |                                                                                                                                        
tACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTcc                                                                                                                                       >  1:157632/1‑140 (MQ=255)
tACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTcc                                                                                                                                       >  2:123325/1‑140 (MQ=255)
 aCCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTccc                                                                                                                                      <  1:1104339/140‑1 (MQ=255)
      cATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTAtt                                                                                                                                 >  1:1012985/1‑140 (MQ=255)
       aTCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTg                                                                                                                                >  1:537437/1‑140 (MQ=255)
            gCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGc                                                                                                                           <  2:157632/140‑1 (MQ=255)
            gCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGc                                                                                                                           <  1:324749/140‑1 (MQ=255)
                tcATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTgg                                                                                                                       >  2:12742/1‑140 (MQ=255)
                 cATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGgt                                                                                                                      >  1:1033794/1‑140 (MQ=255)
                  aTACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGgtg                                                                                                                     >  1:1428276/1‑140 (MQ=255)
                       cTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGc                                                                                                                >  1:593352/1‑140 (MQ=255)
                          aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTgg                                                                                                             >  2:106649/1‑140 (MQ=255)
                          aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCCGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTg                                                                                                              <  2:1433567/139‑1 (MQ=255)
                          aTTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCCGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTg                                                                                                              >  1:1433567/1‑139 (MQ=255)
                            tGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCg                                                                                                           <  1:106649/140‑1 (MQ=255)
                            tGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCg                                                                                                           <  1:12742/140‑1 (MQ=255)
                            tGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCg                                                                                                           <  1:784472/140‑1 (MQ=255)
                             ggCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCgg                                                                                                          <  1:439991/140‑1 (MQ=255)
                               cgcgATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGt                                                                                                         >  2:246726/1‑139 (MQ=255)
                               cgcgATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGt                                                                                                         <  1:246726/139‑1 (MQ=255)
                               cgcgATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTc                                                                                                        <  1:1057899/140‑1 (MQ=255)
                                              tGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTtgctgc                                                                                         >  2:110547/1‑140 (MQ=255)
                                                gTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGt                                                                                       <  2:1479763/140‑1 (MQ=255)
                                                     aCCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGt                                                                                       >  1:884177/1‑135 (MQ=255)
                                                     aCCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGt                                                                                       <  2:884177/135‑1 (MQ=255)
                                                        aaTGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTg                                                                               >  2:1171743/1‑140 (MQ=255)
                                                          tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTTCTGCGTCGCTACTggg                                                                             >  2:316263/1‑139 (MQ=255)
                                                          tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGa                                                                             >  1:433328/1‑140 (MQ=255)
                                                          tGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGa                                                                             >  2:137248/1‑140 (MQ=255)
                                                           gCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAg                                                                            <  2:478467/140‑1 (MQ=255)
                                                           gCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAg                                                                            <  1:185071/140‑1 (MQ=255)
                                                           gCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAg                                                                            <  2:433328/140‑1 (MQ=255)
                                                           gCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAg                                                                            <  2:537437/140‑1 (MQ=255)
                                                           gCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGCATTAAGTATTTTCGGTTCTCTGATCGGCTCCTTCATTGTCCGCCTGGTGTTTGCCGGAGGTCTGATTTTCTGCCTGCGTCGCTAGTgggg                                                                            <  2:423862/140‑3 (MQ=255)
                                                             gTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCg                                                                          <  2:437819/140‑1 (MQ=255)
                                                               tgatgaCCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGc                                                                        >  1:230653/1‑140 (MQ=255)
                                                                    aCCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCaa                                                                   >  1:1479541/1‑140 (MQ=255)
                                                                      cGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGa                                                                 >  1:398119/1‑140 (MQ=255)
                                                                          aCGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAAcg                                                             <  1:1171743/140‑1 (MQ=255)
                                                                             tCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACgcgc                                                          <  2:1516568/140‑1 (MQ=255)
                                                                             tCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACgcgc                                                          <  1:1028534/140‑1 (MQ=255)
                                                                             tCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACgcgc                                                          <  2:1033794/140‑1 (MQ=255)
                                                                             tCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACgcgc                                                          <  1:973930/140‑1 (MQ=255)
                                                                                   gtggATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTAt                                                    >  1:1368801/1‑140 (MQ=255)
                                                                                        tGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAAcg                                                             <  2:718223/126‑1 (MQ=255)
                                                                                        tGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAAcg                                                             >  1:718223/1‑126 (MQ=255)
                                                                                         gCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCt                                              <  2:1428276/140‑1 (MQ=255)
                                                                                         gCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCt                                              <  1:758741/140‑1 (MQ=255)
                                                                                           aCTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGa                                            >  2:1348213/1‑140 (MQ=255)
                                                                                            cTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGAc                                           <  1:137248/140‑1 (MQ=255)
                                                                                             tCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCCGCGTCGCCACTGGAGAGGCACCAAGAAACGCGGCCGTATCCACCCGAcc                                          >  2:758741/1‑140 (MQ=255)
                                                                                               aaTATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGAcccc                                        >  2:229095/1‑140 (MQ=255)
                                                                                                atatCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCa                                       >  2:1117273/1‑140 (MQ=255)
                                                                                                   tCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCg                                    >  2:813823/1‑140 (MQ=255)
                                                                                                       gAAAGTATTAAGTATTTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGGGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCCGCGGAgg                                >  2:387374/1‑139 (MQ=255)
                                                                                                         aagtattaagtattTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGt                              >  1:1479720/1‑140 (MQ=255)
                                                                                                             attaagtattTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGaaa                          >  2:663395/1‑140 (MQ=255)
                                                                                                             attaagtattTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGaaa                          >  2:1475058/1‑140 (MQ=255)
                                                                                                              ttaagtattTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTAGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGaaaa                         >  1:1232673/1‑140 (MQ=255)
                                                                                                               taagtattTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAaaag                        >  1:1154239/1‑140 (MQ=255)
                                                                                                                  gtattTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCGCCTGACCCCAGCGGAGCGTGAaaagaaa                     <  2:1363103/140‑1 (MQ=255)
                                                                                                                  gtattTTCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAaaagaaa                     <  1:1117273/140‑1 (MQ=255)
                                                                                                                        tCGGTTCTCTGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGc               >  1:1160574/1‑140 (MQ=255)
                                                                                                                               tctGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGaaaa        >  2:391981/1‑140 (MQ=255)
                                                                                                                               tctGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGaaaa        >  2:855467/1‑140 (MQ=255)
                                                                                                                               tctGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAAACGGCAAGaaaa        >  2:1517192/1‑140 (MQ=255)
                                                                                                                                ctGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGaaaaa       <  2:593352/140‑1 (MQ=255)
                                                                                                                                 tGATCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAg      >  2:832868/1‑140 (MQ=255)
                                                                                                                                   aTCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCTGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAgcc    <  1:466503/140‑1 (MQ=255)
                                                                                                                                    tCGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAgccg   <  2:1012985/140‑1 (MQ=255)
                                                                                                                                     cGTCTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAgccgc  <  1:123325/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACCGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGC  >  NC_000913/3637988‑3638260

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