Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,760,584 T→G 5.2% L13R (CTA→CGA)  yfjJ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,760,5840TG 5.2% 568.0 / 17.8 154L13R (CTA→CGA) yfjJCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (73/73);  new base G (4/4);  total (77/77)
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.21e-01

Only 100 of 156 total aligned reads displayed.

GGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGC  >  NC_000913/2760455‑2760728
                                                                                                                                 |                                                                                                                                                
ggATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                <  4:687070/148‑1 (MQ=255)
        gCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTa                                                                                                                        >  2:249081/1‑148 (MQ=255)
         cTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAt                                                                                                                       <  1:249081/148‑1 (MQ=255)
             tCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                <  4:989573/135‑1 (MQ=255)
              cATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatg                                                                                                                  >  4:984912/1‑148 (MQ=255)
                           aaTATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCg                                                                                                                                            >  3:1023609/1‑109 (MQ=255)
                           aaTATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCg                                                                                                                                            <  4:1023609/109‑1 (MQ=255)
                             tatataGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCt                                                                                                                                                  >  2:301364/1‑101 (MQ=255)
                                   ggTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTc                                                                                                              <  2:10253/131‑1 (MQ=255)
                                       cATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAt                                                                                                                       <  4:650999/118‑1 (MQ=255)
                                       cATAGACACAGCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAt                                                                                                                       >  3:650999/1‑118 (MQ=255)
                                         tAGACACATCAGAACGATCACTTACTTCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtg                                                                                                                     >  2:488851/1‑118 (MQ=255)
                                         tAGACACATCAGAACGATCACTTACTTCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtg                                                                                                                     <  1:488851/118‑1 (MQ=255)
                                         tAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGaa                                                                                                                                    <  3:99170/103‑1 (MQ=255)
                                          aGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAAt                                                                                      >  3:188116/1‑148 (MQ=255)
                                            acacATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCaa                                                                                       >  2:733028/1‑145 (MQ=255)
                                             cacaTCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGAt                                                                                   >  2:1022357/1‑148 (MQ=255)
                                                 tCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGaa                                                                                                                                    >  3:482028/1‑95 (MQ=255)
                                                  cAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                    <  4:649070/110‑1 (MQ=255)
                                                  cAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                    >  3:649070/1‑110 (MQ=255)
                                                  cAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTAc                                                                              >  1:625202/1‑148 (MQ=255)
                                                    gAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                <  3:1134946/96‑1 (MQ=255)
                                                    gAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                >  4:1134946/1‑96 (MQ=255)
                                                     aaCGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTa                                                                               >  3:272607/1‑144 (MQ=255)
                                                     aaCGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTa                                                                               <  4:272607/144‑1 (MQ=255)
                                                         aTCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGa                                                                       >  3:859746/1‑148 (MQ=255)
                                                          tCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGAt                                                                      >  4:1126316/1‑148 (MQ=255)
                                                           cACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATc                                                                     <  4:188116/148‑1 (MQ=255)
                                                           cACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATc                                                                     <  1:1022357/148‑1 (MQ=255)
                                                             cTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTa                                                                                             <  2:168605/122‑1 (MQ=255)
                                                             cTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTa                                                                                             >  1:168605/1‑122 (MQ=255)
                                                              ttACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                    >  3:989492/1‑98 (MQ=255)
                                                               tACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatg                                                                                                                  >  3:663717/1‑99 (MQ=255)
                                                               tACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATCTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatg                                                                                                                  <  4:663717/99‑1 (MQ=255)
                                                                 cTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATg                                                               <  3:1126316/148‑1 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACt                                                                                          >  2:130701/1‑119 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACt                                                                                          <  1:130701/119‑1 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGAt                                                                      <  3:198713/139‑1 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGAt                                                                      >  4:198713/1‑139 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGAc                                                             >  4:108223/1‑148 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGAc                                                             >  3:306352/1‑148 (MQ=255)
                                                                     aaTATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAtt                                                                                                                      >  3:794450/1‑89 (MQ=255)
                                                                     aaTATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAtt                                                                                                                      <  4:794450/89‑1 (MQ=255)
                                                                     aaTATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGa                                                                                                     <  4:394383/106‑1 (MQ=255)
                                                                     aaTATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTc                                                                                 >  4:1000473/1‑126 (MQ=255)
                                                                      atatCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCaa                                                                                       >  2:592148/1‑119 (MQ=255)
                                                                      atatCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCaa                                                                                       <  1:592148/119‑1 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                    <  1:154634/88‑1 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                    >  1:1111873/1‑88 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                    >  2:154634/1‑88 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                    <  2:1111873/88‑1 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAAt                                                        <  4:859746/148‑1 (MQ=255)
                                                                          cATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                >  4:1080884/1‑74 (MQ=255)
                                                                            ttGGCCAGGATACTATATGTCATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGc                                                                  <  4:749684/134‑1 (MQ=255)
                                                                            ttGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                >  3:859547/1‑72 (MQ=255)
                                                                             tGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCa                                                                                                                                 <  2:315066/70‑1 (MQ=255)
                                                                             tGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCa                                                                                                                                 >  1:315066/1‑70 (MQ=255)
                                                                             tGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATa                                                   <  3:108223/148‑1 (MQ=255)
                                                                               gCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaaa                                                                                                        <  4:977969/93‑1 (MQ=255)
                                                                                  aGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaa                                              >  1:312198/1‑148 (MQ=255)
                                                                                   ggATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGTGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGa                                                              <  4:210580/131‑1 (MQ=255)
                                                                                   ggATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGa                                                              >  3:210580/1‑131 (MQ=255)
                                                                                   ggATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaaa                                             <  4:306352/148‑1 (MQ=255)
                                                                                     aTACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGa                                                     <  3:246840/138‑1 (MQ=255)
                                                                                     aTACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGa                                                     <  4:807536/138‑1 (MQ=255)
                                                                                     aTACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGa                                                     >  4:246840/1‑138 (MQ=255)
                                                                                      tACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAata                                          <  2:625202/148‑1 (MQ=255)
                                                                                            atGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                <  1:767130/56‑1 (MQ=255)
                                                                                            atGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATg                                                               <  3:73815/121‑1 (MQ=255)
                                                                                            atGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGAt                                                    <  4:891514/132‑1 (MQ=255)
                                                                                            atGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAatat                                         <  3:911873/143‑1 (MQ=255)
                                                                                                   aCGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtg                                                                                                                     <  3:761571/60‑1 (MQ=255)
                                                                                                   aCGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTTAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGAt                                                    >  1:1072709/1‑125 (MQ=255)
                                                                                                   aCGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAGCTCAATGATTc                                                                                 >  3:37563/1‑96 (MQ=255)
                                                                                                   aCGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCAtt                             <  3:984912/148‑1 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaa                                                                                                         >  1:701724/1‑69 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaa                                                                                                         <  2:701724/69‑1 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCt                                                                                               <  1:1007407/79‑1 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAAt                                                        >  3:414476/1‑118 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAAt                                                        <  4:414476/118‑1 (MQ=255)
                                                                                                        aGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                >  4:627185/1‑44 (MQ=255)
                                                                                                        aGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                <  3:627185/44‑1 (MQ=255)
                                                                                                           gAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCttt                         <  2:50866/144‑1 (MQ=255)
                                                                                                           gAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCttt                         >  1:50866/1‑144 (MQ=255)
                                                                                                             aTTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaa                                                                                                         <  2:972184/62‑1 (MQ=255)
                                                                                                             aTTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaa                                                                                                         >  1:972184/1‑62 (MQ=255)
                                                                                                             aTTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTg                                                                   >  2:434201/1‑100 (MQ=255)
                                                                                                               tCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCaa                                                         <  3:1017009/108‑1 (MQ=255)
                                                                                                               tCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAg                 >  4:13069/1‑148 (MQ=255)
                                                                                                                       attctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATaa                                   >  2:809202/5‑122 (MQ=255)
                                                                                                                        aTGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCa                                                                                        >  4:603190/1‑68 (MQ=255)
                                                                                                                         tctATCCCGATATAGGGTTAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCaa                                                         >  1:534357/3‑98 (MQ=255)
                                                                                                                         tctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCaa                                                         <  2:534357/96‑1 (MQ=255)
                                                                                                                         tGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGc                                                                  <  1:970591/89‑1 (MQ=255)
                                                                                                                            aTCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCa                                                                 >  3:1081799/1‑87 (MQ=255)
                                                                                                                              cccTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACgt                                                                            >  1:1074730/1‑74 (MQ=255)
                                                                                                                              cccTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACgt                                                                            <  2:1074730/74‑1 (MQ=255)
                                                                                                                              cccTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTg            <  1:802841/138‑1 (MQ=255)
                                                                                                                              cccTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTg            >  2:802841/1‑138 (MQ=255)
                                                                                                                               ccTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATa                                    <  1:45700/113‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                                
GGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGC  >  NC_000913/2760455‑2760728

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